@ Enhanced Elliptic Grid Generation 

Decay parameters that govern grids near boundaries are determined automatically. 
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An enhanced method of elliptic grid 
generation has been invented. Whereas 
prior methods require user input of cer- 
tain grid parameters, this method pro- 
vides for these parameters to be deter- 
mined automatically. 

“Elliptic grid generation” signifies gen- 
eration of generalized curvilinear coordi- 
nate grids through solution of elliptic 
partial differential equations (PDEs). 
Usually, such grids are fitted to bounding 
bodies and used in numerical solution of 
other PDEs like those of fluid flow, heat 
flow, and electromagnetics. Such a grid is 
smooth and has continuous first and sec- 
ond derivatives (and possibly also contin- 
uous higher-order derivatives) , grid lines 
are appropriately stretched or clustered, 
and grid lines are orthogonal or nearly 
so over most of the grid domain. The 
source terms in the grid-generating PDEs 
(hereafter called “defining” PDEs) make 
it possible for the grid to satisfy require- 
ments for clustering and orthogonality 
properties in the vicinity of specific sur- 
faces in three dimensions or in the vicin- 
ity of specific lines in two dimensions. 

The grid parameters in question are 
decay parameters that appear in the 
source terms of the inhomogeneous 
defining PDEs. The decay parameters 
are characteristic lengths in exponen- 
tial-decay factors that express how the in- 
fluences of the boundaries decrease with 
distance from the boundaries. These 
terms govern the rates at which distance 
between adjacent grid lines change with 
distance from nearby boundaries. 

Heretofore, users have arbitrarily 
specified decay parameters. However, 
the characteristic lengths are coupled 
with the strengths of the source terms, 
such that arbitrary specification could 


lead to conflicts among parameter val- 
ues. Moreover, the manual insertion of 
decay parameters is cumbersome for 
static grids and infeasible for dynami- 
cally changing grids. 

In the present method, manual inser- 
tion and user specification of decay pa- 
rameters are neither required nor al- 
lowed. Instead, the decay parameters are 
determined automatically as part of the 
solution of the defining PDEs. Depend- 
ing on the shape of the boundary seg- 
ments and the physical nature of the 
problem to be solved on the grid, the so- 
lution of the defining PDEs may provide 
for rates of decay to vary along and 
among the boundary segments and may 
lend itself to interpretation in terms of 
one or more physical quantities associ- 
ated with the problem. 

The limiting form of the defining 
equations used in this method is partly 
analogous to boundary-value PDEs for 
heat transfer over long, thin fins govern- 
ing convection and conduction of heat 
across the boundary segments and heat 
generated or lost within the volume to be 
enclosed by the grid. Each limiting form 
of the defining equations is deemed to be 
valid near at least one boundary segment. 
Each such equation includes at least two 
independent Cartesian coordinate vari- 
ables and at least one generalized coordi- 
nate as a dependent variable, the integral 
form of which constitutes a boundary 
constraint. Boundary conditions analo- 
gous to temperature and thermal con- 
ductivity prescription can be specified by 
the user. In addition, as an essential ele- 
ment of the method, a selected power of 
at least one heat-transfer coefficient must 
correspond to at least one decay parame- 
ter (which, as stated above, the user does 



A Grid in an Annulus between relatively rotating 
inner and outer circles was generated by the 
present method. 


not specify) near at least one boundary 
segment, and such a heat-transfer coeffi- 
cient evolves according to the boundary 
constraint. 

The figure presents results of applica- 
tion of the method to a two-dimensional 
annular region. This example illustrates 
the clustering of grid points near the 
inner boundary of the annulus. It also il- 
lustrates how the grid evolves automati- 
cally with one or more time-varying 
boundary condition (s) — in this case, 
rotation of the inner boundary relative 
to the outer boundary or vice versa. In a 
case involving time-varying boundary 
conditions, the solution at each time 
step serves as a starting point for the so- 
lution at the next time step. 

This work was done by Upender K. Kaul of 
Ames Research Center. 

Inquiries concerning rights for the commer- 
cial use of this invention should be addressed 
to the Technology Partnerships Division, 
Ames Research Center, (650) 604-2954. 
Refer to ARC-14710-1 


©Automated Knowledge Discovery From Simulators 

Active learning process efficiently explores simulator input space. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A computational method, SimLearn, 
has been devised to facilitate efficient 
knowledge discovery from simulators. 
Simulators are complex computer pro- 
grams used in science and engineering to 
model diverse phenomena such as fluid 


flow, gravitational interactions, coupled 
mechanical systems, and nuclear, chemi- 
cal, and biological processes. SimLearn 
uses active-learning techniques to effi- 
ciently address the “landscape characteri- 
zation problem.” In particular, SimLearn 


tries to determine which regions in “input 
space” lead to a given output from the 
simulator, where “input space” refers to 
an abstraction of all the variables going 
into the simulator, e.g., initial conditions, 
parameters, and interaction equations. 


26 


NASA Tech Briefs, July 2007 


